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We use operator product expansion and QCD light-cone sum rules to estimate the 
amplitude of the decay B — > J/ipK taking into account leading nonfactorizable 
contributions. The result is very similar to the estimate obtained earlier from 
a four-point QCD sum rule. We discuss applications of this method to other 
nonleptonic B decay modes. 



1 Introduction 

Exclusive nonleptonic decays of B and D mesons are complicated processes in- 
fluenced by strong interactions at small and large distances. Theoretically, the 
decay amplitudes are determined from an effective weak Hamiltonian in terms 
of short-distance Wilson coefficients and matrix elements of local four-quark 
operators. While the short-distance effects can be calculated in QCD pertur- 
bation theory using renormalization group methods, it is extremely difficult to 
obtain reliable and accurate estimates for the hadronic matrix elements. Here 
we report on an estimate of the matrix element for the decay B — > J/tpK using 
QCD sum rule techniques. 

The decay amplitude for B — > J/ipK is determined by the matrix element 



{KJ/i; | H w 



b) = ^ 2 v cb v c : 



C2(/«) 



ci(m) 



(KJ/^ | Q 2 (n) | B) 



2 Ci (»)(KJ/tP I 2 (/i) | B) 



(1) 



where Gf is the Fermi constant, V c b and V cs are the relevant CKM matrix 
elements, and Ci(/i) and c 2 (^) are the Wilson coefficients of the four-quark 
operators 



O x = {sTP C ){cT p b), 2 = {cTPc){sT p b), 



(2) 
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respectively, with T p = j p (l — 75). The operator 



with tr(X a X b ) = 25 ab originates from the Fierz rearrangement of 0\. 

In the usual factorization approximation, the matrix element of O2 is split 
into the product 

(KJ/yj I 2 ( M ) I B) = (J/V> I cT'c I Q)(Jf I 5T„& | B) , (4) 
involving simpler matrix elements of quark currents: 

(0 I cT"c I J/V>(p)) = /v-m^e", (5) 

I <T> I B(p + q)) = 2f+(p 2 )q p + (f + (p 2 ) + f~(p 2 ))p p . (6) 

In the above, is the «// 1/> decay constant, e p is the polarization vector, 
and f^(p 2 ) are form factors. The momentum assignment is also indicated. 
The matrix element of O2 vanishes in this approximation because of colour 
conservation. 

The short-distance coefficients ci^in) and the matrix elements entering 
(|l|) are scale-dependent, whereas the decay constants and form factors deter- 
mining the r.h.s. of (|J) are physical scale-independent quantities. Therefore, 
factorization can at best be an approximation valid at a particular scale. As 
shown in ref. u, in the B — ► J/ipK channel factorization does not work at 
fx = Oirrib). The main nonfactorizable contribution comes from the matrix 
element of O2 which is conveniently parametrizedu in the form 

(KJ/ip I 62 (a») I B) = 2f m f B ^ K { i i){e ■ q) . (7) 

Here, the dimensionless quantity /s^j^/x) is the analog of } + {m 2 b ) in (Q). Note 
that the nonfactorizable effects in the matrix element of 0\ are subdominant 
and therefore neglected here. Substitution of ([|) and ([I]) in (|l|) yields 

(KJ/i> \H W \B)= V2G F V cb V* s a^ K Uf+m^e ■ q) (8) 

with 



12 =C2( " ) + T + 2ci( ' i) /TK) 



(9) 



An estimate of ihe matrix element (ffl) using four-point QCD sum rules 
is described in ref. □. In this approach, originally suggested in ref. □ for D 
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decays, gluonic interactions which break factorization are associated with the 
nonperturbative nature of the QCD vacuum as described by quark and gluon 
condensates. The nonfactorizable matrix element (^) turns out to be small in 
comparison to the factorized matrix element ([|). Numerically, one finds 

= -(0.045 to 0.075) (10) 

and f + (m 2 h ) = 0.55. Note that the relevant scale \i in ( |l0| ) is of order of the 
Borel mass M ~ y/m B - m 2 ~ 2.4 GeV. Although the ratio Jb^k(M)/ / + (m|) 
~ 0.1 is small, it has a strong impact in (||) because of the enhancement by the 
large short-distance coefficient c\{M) ~ 1.1 and because of a partial cancella- 
tion in C2(M) +ci(M)/3 ~ 0.09. The resulting prediction for a-i is-jcjomparable 
in size but opposite in sign to the outcome of phenomenological fitsB'ta assuming 
a universal coefficient ci2 (and a\). It should be stressed that in our approach 
the nonfactorizable effects are expected to be nonuniversal. Furthermore, the 
nonfactorizable term in (||) tends to cancel the nonleadipg in 1 /N c factorizable 
term c%/3. This resembles the l/iV c -rule for D decays^ 
An alternative approach was suggested in ref. E 
uct expansion (OPE) and absorbing the nonfactorizable gluon interactions in 
hadronic matrix elements of new effective operators. In ref. El this method was 
applied to the decay mode B° —>■ D + it~. In the heavy quark limit, the new 
operator reduces in this case to the HQET chromomagnetic operator. The ma- 
trix element of the latter between B and D states can be reliably estimated. 
A similar approach can also be applied^ to decay modes where the nonfactor- 
izable amplitudes involve hadronic matrix elements between a B and a light 
meson. Here we describe such an estimate for B — > J/ipK. 



using operator prod- 



2 The correlation function and OPE 

As suggested in ref. i, we replace the J/ ip state in (Q) by the generating current 
C7^c and consider the following correlation function: 

A„(p,q)=i J d^x e^{K(q) \ T{c(x)^c(x) 7 O 2 (0)} \ B(p + q)) 

= (-p 2 q f , + (p-q)p^A(p 2 ), (11) 

where (p + q) 2 — m 2 B and q 2 = m 2 K is put to zero. Because of current conser- 
vation, there is only one invariant amplitude A. For the momentum p of the 
c7 M c current we require 

P 2 <4m^. (12) 
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Inserting a complete set of states with J/ip quantum numbers between the 
operators in (^Tj) one gets, schematically, 

A (p,q)= (°\^ C \ J /^)( KJ M02\B) | (0\c^c\h)(Kh\O 2 \B) ^ 



Substitution of (||) and (0) in ([b^) then leads to the following dispersion relation 
for the invariant amplitude A: 

mi - p 2 ml, - p 1 J s - p^ 

Ami 

The J/i/; ground-state contribution in this relation contains the desired matrix 
element (^). The second term in the r.h.s. represents the contribution of the 
ip' resonance, whereas the dispersion integral takes into account all excited 
charmonium and non-resonant states in J/tp channel located above the open 
charm threshold. Possible subtraction terms are not shown for brevity. They 
will be eliminated later. In the region (|l^), it is possible to expand the product 
of operators entering- (|ll]) in a series of local operators. A similar expansion has 
been applied in ref. LJ in order to estimate the long-distance effect in B — > K*^f 
due to interaction of virtual charmed quarks with soft gluons. As a first step, 
we contract the c-quark fields in the product of the currents cj^c and cT p krC 
and use 

(0|T{ C (x)c(0)}|0) = 



d k ik(x _ y) 



(2tt) 4 1 k 2 



(k 2 - ml) 2 



Ak 2 -m 2 )j x -2k x (fi 
V — 



(fc 2 



(15) 



Here the first term is the free c-quark propagator and the second term in- 
corporates the one-gluon emission in the fixed-point (x — 0) gauge A^(x) — 
1 1 *2x 7 G% ^(0) . The free-quark loop vanishes because of TrX a = 0. The lowest 
order nonvanishing contribution involves a single gluon-field tensor. The cor- 
responding Wilson coefficient can be calculated from the diagrams shown in 
Fig. 1. The result is : 

i f d^xe^Ticixh^ix), cT„y c(0)} = t^ a0 g s G aa0 I c (p 2 ) + ... (16) 



4 



with 



(17) 




Figure 1: Diagrams determining the Wilson coefficients in the OPE (j£ 

The higher-dimensional terms of the expansion denoted by the ellipses in 
( |l6| ) involve operators with derivatives. of the gluon field and/or multiple gluon- 
field operators. As discussed in ref. a, for the physical b and c masses, these 
terms are expected to contribute at the level of 10% or less. We neglect them 
here. 

Substitution of @ in (|ll|) yields 

A(p 2 ) = -9bk(p 2 )I c (p 2 ) , (20) 

where the effective form factor qbk is defined by the matrix element 

t^ af }(K(q) | sr»g s G a "b | B(p + q)) = (p 2 q^ - (p ■ q)p„) g BK (p 2 ) (21) 

with G afj = ^-G aaf3 . Note that the relevant scale for A is the scale of OPE, 
that is /i ~ 2m c . The calculation of the form factor 9bk{p 2 ) requires some 
nonperturbative method. In the next section we will estimate it using light- 
cone sum rules. 

3 Light-cone sum rule for the form factor gsK 



To evaluate the matrix element (|21j) we proceed as in the, calculation of the 
B — > 7T, K form factors in ref. la (see also the reviews ref. Considering the 



5 



correlation function 
Fp(p,q) = it MvaP J d 4 xe ipx (K(q) \ T{s(x)T» g s G al3 (x)b(x), b(0)i l5 d(0)} | 0) 

= (p% ~ (p ■ q)Pv) F( P 2 , (p + q) 2 ) , (22) 

we insert a complete set of intermediate states with B-meson quantum numbers 
between the currents. This yields: 



(K\sT»g s G a Pb\B)(B\bi l5 d\0) 

m B-(p + q) 2 



(K | sVgsG^b | h B )(h B I bi-y 5 d | 0) \ 
Y B m 2 hB -{p + q) 2 J ' 

where the ground state contribution contains the matrix element of interest. 
Using © and 

m b (0 | bi l5 d | B) = m%f B , (24) 

where f B is the S meson decay constant, and representing the sum over excited 
and continuum states by a dispersion integral, one finds 

oo 

F(P 2 , (P + q?) = fl BgB ( Ki f\ 2) + ! ■ (25) 

m b (m z B - (p + q) z ) J s - (p + q) 2 

s h 

In the momentum region [p + q) 2 <C m 2 and p 2 <C 4m;?, the amplitude F can 
also be calculated with the help of OPE. The result can be brought in the form 
of a dispersion integral: 

oo 

F Q cn(p 2 ,(P + q) 2 ) = - Ids Im W*y . (26) 
7T J s-(p + qy 



Furthermore, the spectral density of the higher states in ([25]) is approximated 
by 

p hs (p 2 ,s)e(s- S %) = -lmF QCD (p 2 ,s)e(s~ S B ) , (27) 



as suggested by quark- hadron duality. With (27) it is then straightforward to 
subtract the contribution of the excited and continuum states in the equation 



given by ( p5| ) and (|2^) and solve this equation for gsK- After performing the 
Borel transformation in (p + q) 2 , the solution reads: 



9bk(p 2 ) = 



m b 



nfBm 2 B 



IirFq C d(p 2 ,s) exp 



M 2 



ds 



(28) 



In contrast to conventional sum rules based on the short-distance OPE 
in terms of local operators, here an expansion in terms of nonlocal operators 
near the light-cone, x 2 ~ 0, is used to calculate the invariant amplitude Fqcd- 
Contraction of the 6-quark fields in the correlation function (|||) and use of the 
free 6-quark propagator leads to the diagram Fig. 2a. In this approximation, 
Fqcd is given by a set of bilocal matrix elements (K\s(x)r a g s G a i3(x)d(0)\0) 
where T a denotes a combination of Dirac matrices. These matrix elements are 
parametrized in terms of three-particle light-cone wave functions of different 
twist: 

(K\s(x)a^ 5 g s G a p{x)d{Q)\Q) = if 3 i<[(q 

-(q a q,g f ^ - q q„g a »)] I Va % ^ K ( ai )e lqx ^ +a ^ , (29) 



(K\s(x) ltil5 g s G af3 (x)d(0)\0) = f. 



K 



9/3 g ai > 



qx 



q a gap. 



+f K ^(q a Xf3 - qpx a ) / Va inK {a t y qx ^ +a ^ , (>!()! 



qx 



P a ^±K(ai)e i9K(ai+a3) 



{K\s(x)^g s G af 3(x)d{0)\Q) =if K 



qp g ail , - 



qx 



( xpqA 



+ifK^(q a x ~ q x a ) j Va, !p\\ K (a i )e iqx ^ +a ^ (31) 

with T>a.i = da\da2da 3 8{\ — ot\ — ct2 — a 3 ). The wave function (paxi&i) — 
<P3K (cci, a.2, a 3 ) has twist 3, while <p±k, <P\\k, <P±k and tpujc are all of twist 
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4. In the following, we consider the S f C/(3)-symmetry limit where fax = /37r> 
^3if = <P3tt, <P±k = fi-ni etc. n The explicit expressions for the pion wave 
functions can be found in ref. u. Higher Fock-space components of the K 
meson wave function are neglected as well as perturbative corrections of the 
kind shown in Fig. 2b. Omitting further details of the calculation we present 
the final result: 

iV, (p + q) 2 ) = 2f 3K (lP) I Vat 2 ( Z K ( (a i VY2 

J rn%-{p+ [ai + a 3 )qr 

1 s U 

+2f K m b / — — - — ■ / da 3 (3(pxK(oti) ~ tp\\K(ai)) ai =u-a 3 

J m{- (p + uq) 2 \ J a 2 =i-u 





1 / ml -p 
u \ml — (p + uq) 



U V 

2 „2 



~ 1 ) l ^ ( <f±K< " a ^ + ^WK^i^a^v-aJ -(32) 





Substituting the above result in (ES) yields the light-cone sum rule: 



9bk(p 2 ) = -j-^r / — exp 



l 

2 r du I m% m 2 —p 2 {\ — u) 



f B m% J u ^ \M 2 uM 2 

A 



hn - V3K(oti) + f„ (3^jj r (a i ) - ^if(a«)) L 1=u _ a3 

/ ct 2 =l-« 



+ u("aM^ _1 ) / ^ (^ XiC ( ai ) + ^ll-K( a »))ai=i'-a3 r ' ( 33 ) 

2 J 

where A = to 2 — p 2 /(sq — p 2 ). In the above estimate the scale of the form 
factor gBK is set by the Borel mass M ~ \J m \ — m \- 

4 Estimate of the nonfactorizable amplitude in B — > J/ipK 

From j2C| ) and (|3^) we get an estimate for the invariant amplitude A(p 2 ) which 
is related to Jb-^k by (|l4|). In order to suppress the contribution of excited 
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Figure 2: Diagrammatic representation of the correlation funciton \2ty. Solid lines represent 
quarks, dashed line gluons, wavy lines external currents, and ovals light-cone wave functions 
of the kaon. 



charmonium and nonresonant states in the latter dispersion relation and to 
eliminate subtraction terms we take the n-th derivative of ( |l4| ) at p 2 = 0: 



JB^K + "7- 

U J v^v 



2»+2 



fBip'K 



+ 



„2n+2 

!V 



p h (s) ds m 



2n+2 

1> 



dp 2 



9bk{p 2 )Ic{p 2 



p 2 =0 



(34) 



Since the influence of higher-dimensional operators neglected in (|lq ) increases 
with n we restrict the numerical analysis with n < 8. For the charmed quark 
mass we take m c = 1.3 GeV and for decay constants the values — 405 MeV, 
Id'/fip — 0.7 extracted from the measured leptonic widths of the J/ip and ip 1 
i3. The values of parameters appearing in the light-cone sum rule for gsK are 
given in ref. u. 

The contribution from the ip' to the sum rule (Q) is suppressed relative to 
the J/i[> contribution roughly by a factor 10 at n — 3 and a factor 50 at n — 8. 



It is therefore reasonable to neglect the integral over still heavier states in (34) 
altogether. From the remaining sum rule one can get a first rough estimate of 
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fB^pK by taking n > 6 and dropping the ?p' term. This yields 

/b^kO) = -(0.04, 0.05, 0.07) at n = 6,7,8, (35) 

where /x ~ 2m c = 2.6 GeV which accidentally coincides with the scale M ~ 
^/ — = 2.4 GeV of the form factor qbk used to get (|3|). For a more 
elaborated estimate one may fit the amplitudes /b^k and /b^'k to the mo- 
ments ( gjj ) for n = 2 -j- 8. This gives fBipiiip) — —0.06 and fBip'xip) = +0.3. 

5 Discussion 

QCD sum rule techniques together with light-cone wave functions provide 
new ways to go beyond factorization in exclusive nonleptonic decays of heavy 
mesons. Within this approach, we have estimated the nonfactorizable ampli- 
tude (0) for the decay B — ► J/tpK . The result lEE]) agrees with an earlier 
estimate ( |Io|) obtained from a four-point sum ruleo" Furthermore, it is inter- 
esting to note that the nonfactorizable amplitude Jb^'k for B — > ip'K, differs 
from fBipK both in sign and magnitude. This finding underlines the general 
expectation of nonuniversality of the coefficients a 2 in the sum rule approach. 

Although it is encouraging that both methods give similar results, there are 
still unknown uncertainties. They arise from several sources: neglect of higher- 
dimensional operators, unknown perturbative corrections, a crude model for 
the hadronic spectral density of higher states in the J ftp channel, and numerical 
uncertainties on the parameters and wave functions. The uncertainty on /b^'k 
is obviously larger than the one on fBipK, because the former has a much 
smaller coefficient in ( |3^ ) than the latter. 

Our final comment concerns the original applications in refs. EjlI of this 
method to B — > Dir . Let us consider the B a — ► D°ir° mode. The nonfactor- 
izable amplitude 

(D a Tr° | 6' 2 | B°) (36) 

with 

0' 2 = (cT"^-u)(dT p ^-b) (37) 

has been estimated in ref. I from the correlation function (jr° |T{ 147^75(3, O' 2 }\B ), 
w7ai75C being the generating current of D . Contraction of all quark fields in 
the operator product 

T{S(x) 7 p75c(x),c(0)r"y U (0)} (38) 
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leads to the gluon- field operator, as shown in section 2. However, the OPE of 
( |38| ) also contains the operator u^-T r u resulting from contraction of the heavy- 
quark fields alone. This operator contributes to the sum rule for ( |36|) through 
matrix elements of effective four-quark operators (ir° | (u^-T r u)(dT p ^-^l\ B). 
These contributions which have not been taken into account in refs. Q'Q may 
significantly influence the results for the nonfactorizable amplitudes in B — > 
Dir. 
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